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Abstract 

We present limits on ultra-high energy (UHE; E v > 10 15 eV) neutrino fluxes from gamma-ray 
bursts (GRBs), based on recently presented data, limits, and simulations from the RICE experi- 
ment. We use data from five recorded transients with sufficient photon spectral shape and redshift 
information to derive an expected neutrino flux, assuming that the observed photons are linked 
to neutrino production through pion decay via the well-known "Waxman-Bahcall" prescription. 
Knowing the declination of the observed burst, as well as the RICE sensitivity as a function of 
polar angle and the previously published non-observation of any neutrino events allows an estimate 
of the sensitivity to a given neutrino flux. Although several orders of magnitude weaker than the 
expected fluxes, our GRB neutrino flux limits are nevertheless the first in the PeV-EeV energy 
regime. For completeness, we also provide a listing of other bursts, recorded at times when the 
RICE experiment was active, but requiring some assumptions regarding luminosity and redshift to 
permit estimates of the neutrino flux. 
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I. INTRODUCTION 



A. Gamma-ray Bursts 

Gamma-ray bursts (GRBs) are the most energetic explosions in the Universe, releasing 
upwards of ~ 10 51 ergs of radiation energy in 7-rays (~ 0.1-1 MeV) within a few seconds (see 
Refs. [1, 2] for reviews on GRB observations). The Burst and Transient Source Experiment 
(BATSE) on board the Compton Gamma-Ray Observatory (CGRO) satellite detected bursts 
at the rate of about 300 per year, in the period 1991-2001, isotropically distributed over the 
sky. The observed duration of the bursts falls into two categories, long bursts lasting for 
tens of seconds and short bursts lasting for less than 2 seconds. The BeppoSAX, High 
Energy Transient Explorer (HETE-2) and SWIFT satellites are observing GRBs currently 
on a routine basis. However, the redshift measurement, and hence the total energy estimate, 
is possible only for a small fraction of the bursts as it requires identification of the GRB host 
galaxies in optical wave bands. 

The unsurpassed 7-ray luminosity naturally make GRBs candidate sources of the observed 
ultra-high energy cosmic-rays [3, 4], many with > 10 18 eV. High energy [E v > 100 GeV) 
neutrinos may then be created from the cosmic-rays, as in beam-dump experiments, while 
propagating in the cosmic photon backgrounds [5-7] and/or at the particle acceleration sites 
of the sources [8-10]. While diffuse flux estimates are a useful method for calculating the total 
neutrino energy injection into the universe by GRBs, individual source detection is crucial for 
understanding the source production processes, and in particular, for verifying whether point 
sources are the sources of UHE cosmic rays. For transient sources such as GRBs, individual 
luminous and bright bursts may dominate the observed diffuse emission as well. A number 
of authors have calculated expected neutrino production from individual long bursts (short 
bursts are typically lower in energy output and hence less likely to be detected individually in 
neutrino telescopes[ll-16]). When folded in with a given experimental acceptance, such flux 
predictions can then be translated into expected sensitivities or, in the case of no detection, 
upper limits. 



B. Radiowave Neutrino Detection 

The RICE experiment has goals similar to the larger AMANDA and ICECUBE experi- 
ments - both seek to measure UHE neutrinos by detection of Cherenkov radiation produced 
by neutrino-nucleon interactions. Whereas AMANDA/ICECUBE is optimized for detection 
of penetrating muons resulting from + N — > fx + N', RICE is designed to detect com- 
pact electromagnetic and hadronic cascades initiated by charged or neutral current neutrino 
interactions in a dense medium: vi(/Vi) + N — > ^ + N' or vi{/v{) + N — > vi{/Vi) + N' 
respectively, with I — e, //, r. As the cascade evolves, a charge excess develops as atomic 
electrons in the target medium are swept into the forward-moving shower via Compton scat- 
tering and positrons are depleted via annihilation, resulting in a net charge on the shower 
front of Qtot ~ E s e/A; E s is the shower energy in GeV [17-24]. Such cascades produce 
broadband Cherenkov radiation - for Ae^i^ 3> r Molicre , the emitting region approximates 
a point charge of magnitude Q tot and therefore emits fully coherently; fortuitously, the field 
attenuation length at such wavelengths in cold polar ice (~ 1 GHz) exceeds 1 km [25]. The 
quadratic rise of the Cherenkov power with neutrino energy, combined with the long attenu- 
ation length lead to the expectation that the effective volume of radio-based detectors begin 
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to exceed those of photomultiplier-tube based detectors at a cross-over energy of ~1 PeV 
[26], leading to neutrino detector arrays consisting of tens of buried radio antennas. The 
RICE radio antenna detection concept of the Cherenkov radiation produced by an in-ice 
weak current event is depicted in Figure 1. 

In this paper we determine model-dependent RICE neutrino flux upper limits for five 
GRBs, namely 020813, 020124, 050603, 050730 and 050908 for which we have extensive ob- 
servational data including redshifts. Using calculations published elsewhere of the expected 
radio-frequency signal strength due to an electromagnetic shower [18, 19], the RICE hard- 
ware, reconstruction software and simulation [27-29], and a neutrino analysis based on data 
taken 1999-2005 [29], we can translate an expected neutrino flux from the bursts recorded 
while RICE was active into an expectation for the RICE neutrino yield. These expectations 
are valid in an energy regime somewhat higher than those typical of, e.g., previous AMANDA 
gamma-ray burst studies [12]. Our calculation of expected neutrino fluxes from these GRBs 
is the most novel aspect of this analysis, and is presented in greater detail herein. 

II. HARDWARE 

The RICE experiment presently consists of a 20-channel array of dipole radio receivers, 
scattered within a 200 mx200 mx200 m cube, at 100-300 m depths. The signal from each 
antenna is boosted by a 36-dB in-ice amplifier, then carried by coaxial cable to the surface 
observatory, where the signal is filtered (suppressing noise below 200 MHz), re-amplified 
(either 52- or 60-dB gain), and split - one copy is fed into a CAMAC crate to form the event 
trigger; the other signal copy is routed into one channel of an HP54542 digital oscilloscope. 
Four antennas which register "hits" within a 1.25 microsecond window define an event trig- 
ger. "Hits" are defined as a voltage excursion greater in magnitude than a discriminator 
threshold, determined at the beginning of each run as that threshold which gives approxi- 
mately 90% livetime. Deadtime is incurred by a full event trigger readout (approximately 
10 s of deadtime), the fast software rejection of anthropogenic noise (approximately 10 ms 
of deadtime per calculation), or the fast hardware rejection of anthropogenic noise (1.2/xsec 
per event). Short- duration pulses broadcast from under-ice transmitters provide the primary 
calibration signals and are used to verify vertex reconstruction techniques. 

Given the known experimental circuit gains and losses [27], the effective volume V e ff is 
calculated as a function of incident E u , as an exposure average of the detector configurations. 
The most important variable is the global discriminator threshold, which is adjusted to 
maintain an acceptable trigger rate under conditions of varying environmental noise. 

Additional details on the RICE shower search procedures are presented elsewhere [29]. 
The most recent search has resulted in zero neutrino candidates, yielding competitive limits 
on the incident ultra-high energy neutrino flux at energies above 10 17 eV. 

III. GRB MODEL AND NEUTRINO FLUX CALCULATION 

In the most widely accepted GRB model, known as the fireball shock model, the prompt 
7-rays are produced by collisions of plasma material moving relativistically (with a bulk 
Lorentz factor T > 100) along a jet (internal shocks), i.e. a fireball. Late time collisions of 
jetted material with an external medium (external shocks) produce X-ray, UV and optical 
radiation, collectively known as GRB afterglow (see Ref. [30, 31] for reviews). A signifi- 
cant number of nucleons are expected to be present in the jet along with leptons (electrons 
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FIG. 1: RICE concept. Shown is an incident neutrino interacting ^100m below the surface and 
illuminating the RICE radio antenna array, indicated by spheres and drawn to scale. Horizontal 
extent is approximately 4 00 m. 

and positrons). The probable mechanism(s) responsible for the observed photons is/are 
synchrotron radiation or/and inverse Compton scattering by high energy electrons. These 
electrons are stochastically accelerated by relativistic shocks via the well-known Fermi mech- 
anism in the tangled local magnetic field resulting in a power-law energy distribution. 

Protons are expected to co-accelerate with electrons in the shocks. High energy non- 
thermal neutrinos are expected to be produced by photonuclear (jry) interactions of protons 
with observed prompt 7-rays [8, 9] and afterglow photons [10] in these optically thin (for 
Thomson scattering) environments. In this analysis, we put limits on the Waxman-Bahcall 
GRB neutrino burst and neutrino afterglow flux models, which are based on the standard 
fireball-shock model of GRBs. We do not consider in the present work other neutrino flux 
models, e.g., (i) precursors from GRB jets buried under the stellar surface [32, 33]; (ii) GRB 
blast-wave models [34] in which both prompt and afterglow emission arise from the external 
shocks; (iii) afterglow from GRBs in dense stellar winds [35]; (iv) the supranova model [36- 
38] in which a supernova type remnant shell from the progenitor star is ejected prior to the 
GRB event. Precursor neutrinos [32, 33] below the RICE threshold [E v > 10 PeV) may 
be detected by neutrino telescopes (e.g. IceCube) without any detectable electromagnetic 
counter-part from their optically thick environment. 

Below, we detail the prescription we have employed in estimating the expected neutrino 
fluxes from a GRB with a well-measured 7-ray flux, following closely the calculations of 
Waxman and Bahcall ([8-10], see also [39, 40]). 

A. GRB observed and model parameters 

The isotropic-equivalent 7-ray bolometric energy and luminosity of a GRB at redshift 
z, given the measured bolometric 7-ray fluence S 1 (typically measured over one decade in 
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TABLE I: Observed and estimated GRB parameters 



GRB 


S, 




*90 


z 


d L 




a (3 


E ■ 




RICE 




(ergs/cm 


2 ) 


(s) 




(Gpc) 


(keV) 




(ergs) 


(ergs/s) 


< Live T > 


050908 


5.1 x 10" 


-7 


20 


3.35 


29.02 


21* 


- 1.93 


1.2 x 10 52 


2.6 x 10 51 


0.93 


050730 


4.4 x 10- 


-6 


155 


3.97 


35.52 


71* 




1.3 x 10 53 


4.3 x 10 51 


0.75 


050603 


3.4 x 10" 


-5 


10 


2.82 


23.59 


289 


0.79 2.15 


5.9 x 10 53 


2.3 x 10 53 


0.92 


020124A 


6.1 x 10- 


-6 


46 


3.20 


27.47 


87 




1.3 x 10 53 


1.2 x 10 52 


0.83 


020813A 


8.4 x 10- 


-5 


89 


1.25 


8.71 


142 


- 2.15 


3.4 x 10 53 


8.6 x 10 51 


0.98 



energy range) and duration t 90 , may be calculated as 

£ 7 ,iso = 47rd 2 S' 7 /(l + z) ; L 7)iso = A7rd 2 L S y /t 90 , (1) 

where is the luminosity distance. The actual jetted energy is (4ir/Q)E ltiso , where fl is 
the jet opening solid angle (typically ~ 0.1 radian). The 7-ray energy, arising from shock- 
accelerated electron synchrotron radiation and/or inverse Compton scattering, is assumed 
to be a fraction e e < 1 of the isotropic-equivalent kinetic energy, i?kin,iso, of the relativis- 
ts plasma material in the standard fireball-shock model. The shock magnetic field shares 
another fraction e B < 1 of -Ekin,iso- 

The luminosity distances may be calculated using the conversion calculator provided by 
the NASA Extra galactic Database (NED [41]) for a ACDM fiat cosmology with H = 70 km 
s~ l Mpc -1 , Vt m = 0.3 and Q A = 0.7 (1 Gpc = 3.1 x 10 27 cm). We have tabulated different 
observed quantities in Table I from 5 GRBs. The calculated values of -E 7i i S o and L 7i i SO are 
also tabulated for later use. The average fractional RICE livetime over the duration of the 
burst is also indicated in the last column. 

The observed 7-ray spectrum from a GRB may be approximately fitted with a broken 
power-law (Band fit [42]) as 

^ocK ; 67 < e ^ b (2) 

G?6 7 I ^ 7 ; £7 ^ ^7,6' 

The spectral indices and break energy e 7) & from observations (where available) are listed in 
Table I. Values with asterisks have been estimated using the phenomenological Ghirlanda 
relationship [43] as 

^T6 = — f-^-V' 68 keV - ( 3 ) 

7 ' 1 + z V10 53 ergs/ v ; 

Note, however, that the precise value of e 7i b is not necessary to estimate the expected neutrino 
flux in the energy range of our interest, as we calculate in the next subsection. 

B. Burst neutrinos 

The condition for resonant p7 interaction at A + energy is e' p e' y ~ 0.3 GeV 2 in the local 
rest frame of the plasma [We denote primed (unprimed) variables in the plasma rest (lab) 
frame. The obesrved variables are redshift-corrected lab variables]. On average, the proton 
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loses a fraction <x p ^&+>— 0.2 of its energy in this interaction. The secondary A + , produced 
by a £>7 interaction decays as either A + — >• nn + or A + — > jot . [53] The pions further decay 
as 7r + — > — > e + u e u fl u^ 1 or 7r° — >• 77. In what follows, we assume equipartition of 
energy among the 4 final state leptons from 7r + -decay. The observed neutrino break energy 
corresponding to the A + resonance is then 

0.015rf / e 76 X- 1 ^ w 



(1 + z) 2 VGeV, 

where e 7) & is the break energy in the 7-ray spectrum [or the peak energy in the e 7 (diV 7 /de 7 ) 
energy spectrum]. The synchrotron break energy (above which pions lose significant energy 
by synchrotron radiation before decaying to neutrinos) in the neutrino spectrum is 



10 llp /™n-J 




TTT—jl^) GeV, (5) 



where the magnetic field in the GRB fireball is given by 



1/2 / "p \ —3 / X v -1 



B - = = 5 >< 104 (life) (&)" (»3h; «• <«) 

The typically unknown parameters are e e « « 0.1; the Lorentz factor of the GRB fireball 
is Tj 300 from various fits. The radius of the fireball is ~ 2F 2 ct v , which varies over a 
large range depending on Tj and the variability time scale (typically, 0.001 s < t v < 1 s). 

To calculate the expected neutrino flux from the GRB fireball, we assume a shock- 
accelerated proton luminosity in the fireball L Piiso ~ L 7)iso with energy distributed as a 
power-law: dN p /de p oc e p 2 up to the maximum energy < 10 20 eV. The total proton energy, 
equal per decade in energy, is then ~ 10-E 7) i SO the 7-ray energy; roughly consistent with the 
assumed equipartition fraction e e ~ 0.1 of £kin,iso. Due to interactions of these protons with 
the observed 7-rays the corresponding neutrino energy flux from the GRB prompt phase, 
following Eqs. (2), (4) and (5), is given by: 



a-;; = \ f -^- x < 
too 



~" " 2 4 to 



v,b\ 



[^vl^v,b\) a 1 ; £v,sb > > Ci/,61 (7) 
{€v,sb / €v,bl) {^u/^u,sb) i £j/ ^ ^v,sbi 

for each neutrino flavor z/ M , and z/ e at the source. The factors 1/2 and 1/4 arise due to the 
WB assumption of equal decay probabilities: A + — > 7r + /7r° and energy equipartion among 
the 7r + decay products, respectively. Here f n is the conversion efficiency of shock accelerated 
protons to pions at the A + resonance in the GRB fireball. For optically thin sources, such 
as GRBs, f n < 0.2. The detailed calculation of f n , which involves many unknown variables, 
is not relevant here. Moreover, the internal shocks (and thus 7-ray emission) occur over a 
wide range of r\ which in turn affect f n . We assume a fixed value f n = 0.2 for our calculation 
which corresponds to a jry — > A + optical depth of the order of unity. The prefactor in the 
expression for flux [Eq. (7)] can be written as 



A j,^ = 5 -7 \ U2HL\ GeV cm" 2 s" 1 , (8) 

8 t 90 V - 2 / V 10 ergs/cm 2 J VlO sj ' V ; 

and the 7-ray spectral indices are a — 1, (5 — 2 to good approximation (1 erg = 624 GeV). 
After oscillations in vacuum, the neutrino fluxes for different flavors are <£v e +p e = = 
&u T +i> T = each with a duration of tburst = £90- We have listed values for these parameters 
in Table II for 5 bursts. 
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TABLE II: Burst neutrino flux parameters 



GRB A (GeV cm' 2 s" 1 ) e„, a6 (GeV) e uM (GeV) t hmst (s) 

050908 4.0 x 10" 7 9.5 x 10 7 3.4 x 10 6 20 

050730 4.4 x 10" 7 6.4 x 10 7 7.7 x 10 5 155 

050603 5.3 x 10" 5 1.2 x 10 7 3.2 x 10 5 10 

020124A 2.1 x 10" 6 4.6 x 10 7 8.8 x 10 5 46 

020813A 1.5 x 10" 5 1.0 x 10 s 1.9 x 10 6 89 



C. Afterglow neutrinos 

The GRB afterglow, arising from external forward and reverse shocks, takes place at a 
radius r e = (3£ , kin,iso/[47rn e x"7.pC 2 r 2 ]) 1 / 3 , which is much larger than r,. Here we assume the 
total jet kinetic energy -Ekin.iso = E lt i so /e e ~ 10-E 7i i SO throughout our calculation and n cx = 1 
atom cm" 3 is a typical interstellar medium (ISM, dominantly neutral hydrogen) density. 
The bulk Lorentz factor of the reverse- shocked plasma shell at the time t = tg when the 
reverse shock has crossed the shell is 



6 4 



17 Ek 



kin,iso 



7rn ex Tn p c 5 tgo 



1/8 — feV'V^rv^r 8 . w 



10 54 ergs J V10 sj Vl cm 



2U/2 



The reverse shock magnetic field in the comoving plasma frame is B' e = (S7re B n ex m p c 
and we set Eb = 0.1 as in the internal shocks. The minimum Lorentz factor of the reverse 
shock-accelerated electrons in the un-shocked plasma frame is 7^ min ~ e e {m p /m e ){Yi/Y e ) and 
the peak electron synchrotron photon energy is e' im — hcT e (3'f' 2 min eB' e )/(2m e c 2 ). However 
the characteristic shock-accelerated electron Lorentz factor, from equating the synchrotron 
cooling time to the dynamic time: td yn = r e/4r 2 c, is r )' ec = 6nm e c/ (<7Th-B e 2 r 4 £dyn) ■ The 
corresponding peak synchrotron photon energy radiated by these characteristic electrons in 
the same B' e field is then 



n r H%eB' e n a ( gj^iso V 25/24 ( r, \ 4/3 ( Y /8 ( n - V 11/24 v 
= hcTe ^^~ OA [io^ vsooJ lioiJ [T^J eV ' (10) 



The number density of these photons in the reverse shock region is approximately given by 
the total synchrotron radiation power (typically at e^ m ) by all electrons scaled at ^ and 
then divided by ^ as 



/ ^ -Ekin,iso Z^ 6 ^ ^ ^ e ( ^7i m 

7,c 8ir 2 hcr 2 \ Tj / m e c 2 m p c 2 \ e' lc 



1/2 



6.8 x 10- (-^Y"" (i) 2/3 f^)' 3/1 f-!=4) 19/12 cm- (11) 
ll0 54 ergs/ V300/ U0 s/ VI cm" 3 / v ; 



10 54 ergs 

The corresponding photon number spectrum is given by 

<^ - < c x { w<°KT ; < < < < c (12) 
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TABLE III: Afterglow neutrino flux parameters 



GRB 


U 


B (GeV/cm 2 /s) 


eu,b2 (GeV) 


e^max (GeV) 


^glow(^) 


050908 


0.03 


6.9 x 10~ 8 


1.5 x 10 8 


1.4 x 10 8 


20 


050730 


0.2 


4.4 x 10~ 7 


2.5 x 10 8 


2.8 x 10 8 


155 


050603 


0.2 


5.3 x 10~ 5 


1.0 x 10 10 


6.0 x 10 8 


10 


020124A 


0.2 


2.1 x 10~ 6 


7.4 x 10 s 


3.3 x 10 8 


46 


020813A 


0.2 


1.5 x 10~ 5 


2.0 x 10 9 


8.4 x 10 8 


89 



Note that, one needs to explicitly calculate the proton to pion conversion efficiency factor 
f n in the afterglow phase (r e > r^) as it may be lower than 0.2, contrary to the prompt 
(internal shocks) phase. We calculate /V for reverse shock-accelerated protons to interact 
with photons of energy e' [Eq. (10)] at the A + resonance as 

F, \ 33 / 24 / + \ -9/8 / „ x 9/8 

^kin.iso \ / l 90 \ / n <™ \ 



/. = mind. ^ + )<^>« 0.2^^ j , (13) 

where Tp 7 ^ A+ ~ cr pi -*&+(j' e /T e )n' ryfi is the optical depth using Eq. (11) and cr pi ^ A + = 
5 x 10~ 28 cm 2 is the py cross-section at the A + resonance. The neutrino break energy for 
the characteristic photon energy [Eq. (10)] is 

e ^ = TTMl&vj GeV - (14) 

The maximum shock-accelerated proton energy can be found by equating the shock- 
acceleration time to the shorter of the synchrotron cooling time (most efficient energy loss 
channel) and the dynamic time. The corresponding maximum neutrino energy is 



e. 



v, max 



<X p _&+> 

4(1 + z) 



(eT e B' e r e ) (15) 



in the afterglow phase. The resulting afterglow neutrino spectrum, following Eqs. (12), (13) 
and (14), is 

e 2 $ s = fnS, x f (£„/£,,,„) ; c < e„, b2 

,62/ j ^,max _ _ ^v,b2i 

for each neutrino flavor and v e at the source. Again, the flavor ratios on Earth, after 

vacuum oscillation, would be $ Ue +u e = = $z/ t +p t = ®i, each with a duration of 

£ g iow = ^90- The parameters /„-, e v ^2 and ej, )inax are listed in Table III for each GRB. The 
pre-factor in the neutrino spectrum B = (/ 7r /8)(S , 7 /tgo), for different GRBs, is also listed in 
Table III. 

The formalism outlined above is schematically presented in Figure 2, which qualitatively 
indicates the relationship between the neutrino and the photon spectra. 



D. GRBs with incomplete data 

We note ~100 GRBs (listed in the Appendix) having insufficient electromagnetic data 
recorded while RICE was live. The break energy (e 7) &) and spectral indices (a and 0) of 
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(a) 




; > log (E v ) 

•v,b 



FIG. 2: Cartoon of the neutrino energy spectrum expected from a GRB having an observed j-ray 
spectrum, (a) Shock accelerated proton spectrum: oc e~ q , and broken power-law fit to the j-ray 
spectrum: oc e~ a for e 7 < e 7i b and e 7 > e^^. (b) Opacity for p^y — > A + resonant interaction 
satisfying the condition e p e 7 = 0.3r 2 GeV 2 . Note that the optical depth changes its shape when 
expressed as a function of e p (lower curve) rather than e 7 (upper curve) to satisfy the resonance 
condition e p e 7 = const, (c) The neutrino spectrum follows the proton spectrum modulo the shape 
of the optical depth. At high energy (not shown here) above the synchrotron break energy [Eq. (5)], 
the neutrino spectrum would steepen by a factor of 2 in the spectral index. 

the 7-ray spectrum may be estimated [e.g., e 7) & from Eq. (3)] and/or sampled from their 
observed distributions. The relative error using the estimated spectral information is small 
and does not largely affect the neutrino flux predictions in the energy range where RICE is 
sensitive. However, the most important piece of information that is lacking is the redshift 
measurement. The luminosity and hence energy of a burst [see Eq. (1)] depends on the 
redshift or its luminosity distance. These in turn affect the model neutrino flux calculation. 
Various methods have been proposed to indirectly measure the redshift of a burst, e.g., 
from the time lag between the hard and soft 7-ray arrival from a burst [44, 45]; using an 
empirical relation between the observed variability and the luminosity [46, 47]; using observed 
correlation between 7-ray total and peak energy [48] ; etc. For now, we have focused on those 
bursts for which reliable redshift data exist. 

IV. SEARCH FOR GAMMA-RAY BURST COINCIDENCES 

Software algorithms used to winnow the RICE data set are presented elsewhere [27-29]. 
Typically, antenna 'hits' must exceed a threshold corresponding to approximately 5.5<r rms , 
where u rms is the rms noise in the recorded waveform for a given antenna channel. In cases 
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Acceptance vs theta and discriminator threshold 



FIG. 3: Relative dependence of V e ff on nadir angle (i.e., ir 7 '2 + declination) , and discriminator 
threshold, based on simulations of the RICE detector, ir /2 corresponds to flux incident along the 
horizon; n corresponds to flux incident from zenith. 

for which the event time is known a priori, one can employ more sophisticated, albeit more 
time-consuming event reconstruction techniques. These include use of a 'matched filter' 
to search for lower-amplitude signals, thereby enhancing the effective volume. Although 
transients such as GRBs lend themselves to such approaches, we have not attempted such 
an analysis and instead rely on the results of previous RICE analyses for sensitivity and 
efficiency estimates. 

We select gamma-ray bursts which were recorded within a ±1000 s overlap time window 
during which the RICE experiment was active [52]. Only GRBs which have declination 
angles less than zero and have therefore been localized to the southern hemisphere are con- 
sidered, given the opacity of the earth to neutrinos at ultra-high energies. The recorded 
RICE experimental discriminator threshold allows a determination of the effective volume 
specifically for the run during which the gamma-ray burst was observed. The record of trig- 
gers, and software surface-noise vetoes over the 2000 second coincidence window allows an 
estimate of the average livetime fraction over that period. Relative to the published RICE 
effective volume and efficiency [29], our ultimate sensitivity is degraded by this calculated 
average livetime fraction. 

Relative to the effective volume integrated over angle [29], appropriate for an expected 
isotropic (i.e., flat in both dN/d(cos8)) and dN/d<f>) diffuse neutrino flux, we must determine 
the acceptance for point sources as a function of declination and azimuth. Our simulations 
indicate that the acceptance is constant (to within ~10%) in azimuth and that the dominant 
angular dependence is on the polar angle of the point source. Since the RICE array is 
shallow, the acceptance for neutrinos incident from the zenith (directly overhead) is nearly 
zero. Maximal acceptance is achieved for neutrinos incident closer to the horizon interacting 
in the 2.5 km of ice below the array which illuminate the array with a relative angle between 
the neutrino momentum vector and the vector to the array close to the Cherenkov angle. 
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FIG. 4: Dependence ofV e ff on nadir angle, integrated over energy. Also shown is the distribution 
of the expected incident flux. The overall efficiency is therefore the properly weighted average of the 
polar angle dependence. 



TABLE IV: RICE limits and ratios to the predicted burst and afterglow fluxes at a given energy 



GRB 


e„ (GeV) 


Burst 
Limit (GeV/cm 2 ) 


Ratio 


e v (GeV) 


Afterglow 
Limit (GeV/cm 2 ) 


Ratio 


050908 


1.0 x 10 8 


20152.48 


2.01 x 10 9 


1.0 x 10 8 


44487.72 


3.21 x 10 10 


050730 


1.0 x 10 8 


24624.48 


6.35 x 10 8 


2.0 x 10 s 


25037.32 


3.04 x 10 8 


050603 


1.0 x 10 8 


8613.60 


8.13 x 10 8 


2.0 x 10 8 


2739.27 


1.68 x 10 8 


020124A 


1.0 x 10 8 


15436.97 


5.44 x 10 8 


2.0 x 10 8 


12874.04 


3.24 x 10 8 


020813A 


1.0 x 10 8 


13975.23 


7.55 x 10 6 


3.0 x 10 s 


2139.85 


7.05 x 10 6 



Figure 3 shows the relative effective volume as a function of the discriminator threshold 
and the nadir angle of the incident neutrino. Given the fact that the angular dependence 
is more important than the energy dependence, Figure 4 shows the angular dependence of 
the effective volume, integrated over energy. For each GRB, this graph provides a correction 
relative to the previously-published V e ff [29] using the tabulated declination angle for each 
analyzed GRB. 

The non-observation of any in-ice shower candidates, coupled with the tabulated RICE 
neutrino-finding efficiency, allows an estimate of the expected event yield for GRB neutrino 
fluxes calculated as above. Our derived upper limits, following the prescription for relating 
a flux model to an upper limit as outlined in the previous RICE neutrino flux study [29], 
are presented in Figure 5. We also list the limit values for each GRB in Table IV and their 
ratios to the predicted flux values at a given neutrino energy. Simple linear scaling then may 
be used to get limit values at another energy. 

An alternative approach to setting limits on the neutrino flux from individual GRBs is to 
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FIG. 5: Upper limit, derived from five indicated GRBs and for SGR 1806-20, assuming a 1:1:1 
isoflavor mix at the detector. Lower set of lines indicate predicted GRB and SGR 1860-20 neutrino 
fluxes as described in the text; higher set of lines indicate RICE upper limits using previously 
published effective volumes and the expected angular dependence of the acceptance. 

use a source stacking method. These methods have been applied in gamma-ray astronomy 
and were discussed recently in the context of neutrinos from various classes of active galactic 
nuclei (AGNs) [49]. The total neutrino flux from all five GRBs in our sample, for the burst 
and afterglow models described, is displayed in Figure 6. We derive upper limits on these 
total fluxes by following the same prescription as for the single burst analysis, using the sum 
of the number of events expected from each individual GRB. 

A class of low luminosity, low energy (dominantly in the X-ray bands) bursts are known 
as X-ray flashes (XRFs). The typical break energy e 7i b ~ 1-10 keV implies that protons 
need to be accelerated to ^ 10 18 eV to produce neutrinos at these sources which is unlikely 
because of their low luminosity. Moreover, the proton to pion conversion efficiency f n would 
be much lower in this case. There were 3 XRFs recorded while RICE was active namely 
GRBs 020903A, 030429A and 030723A (see Appendix) with bolometric fluences 5 7 one to 
two orders of magnitude lower than typical long GRBs listed in Table I. 
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FIG. 6: Upper limit on the total stacked flux from the five GRBs, 020813, 020124, 050603, 050730 
and 050908, assuming a 1:1:1 isoflavor mix at the detector. The lower pair of lines indicate the 
predicted total neutrino flux from the five GRBs for the burst and afterglow models for neutrino 
production; the higher set of lines indicate RICE upper limits for the total neutrino fluxes using 
previously published effective volumes and the expected angular dependence of the acceptance. 

A. Other Transient Source Searches 

Soft Gamma-ray Repeaters (SGRs) represent one class of extreme X-ray pulsars repeat- 
edly emitting ~ 0.1 s bursts of soft 7-rays. SGRs are highly magnetized (~ 10 15 G) neutron 
stars, often called 'magnetars'. Four SGRs have been identified in the Milky Way (1806-20, 
1627-41, 1900+14, and 1801-23) and one has been identified in the Large Magellanic Clouds 
(0526-66). Giant flares from an SGR, hypothesized to originate from a sudden shift of mag- 
netic field, occur with periods of ~ 10 years. In particular, the giant flare of SGR 1806-20 
on 27 December 2004 at 21:30:26 UT was the brightest cosmic transient observed to date 
with a flux ~ 10 ergs s -1 cm -2 lasting for ~ 0.1 s. We note that this flux is ~ 10 4 times 
larger than the most luminous GRBs ever detected. Detected radio afterglow of SGRs imply 
the presence of relativistic jets, similar to GRBs. If high energy neutrinos were emitted 
with comparable energy emitted in 7-rays, the 27 December 2004 giant flare of SGR 1806-20 
might have been detected by operating neutrino telescopes, as predicted by several authors 
[50, 51]. 

We investigated the possibility of a coincidence with the 27 Dec 04 SGR 1806-20 giant 
flare. Although the RICE experiment was active during that time, we did not observe a 
trigger at the time of the SGR flare. We did record two forced triggers ("unbiased" events) 
at 21:25:33 UT and 21:35:33 UT, respectively, which bracket the flare time by almost exactly 
five minutes. A channel-by-channel comparison of the time-domain and frequency- domain 
waveforms for these events reveals no obvious differences between the data recorded before- 
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flare vs. after-flare. Given no detection we then put limits on two model fluxes 




GeV cm-V 1 (Model A) 



(17) 



where Model A corresponds to the neutrino flux originating from pp interactions in Ref. [50] 
and Model B is a generic pj flux, similar to the prediction in Ref. [51], assuming total proton 
energy is equal to the observed 7-ray energy. The maximum neutrino energy in both cases 
is 3.3 x 10 8 GeV [50]. The flux predictions and the limits are plotted in Fig. 5. The limits 
are 191 and 550 times higher than the flux predictions for Models A and B respectively. 

V. SUMMARY 

Selecting GRB events recorded while the RICE experiment was actively taking data, and 
using previously presented results on the RICE analysis procedures and simulations, we have 
estimated the expected RICE sensitivity to such transients. Based on the lack of observed 
coincidences, we set limits on the UHE neutrino flux for several candidates. Although the 
limits presented are well above model predictions, they do, nevertheless, probe possible 
neutrino fluxes from GRBs in an energy regime considerably higher than those previous 
studied. Future improvements in sensitivity will be realized by enlarging the scale of the 
radio detector array, as well as improved individual GRB parameter measurements. 
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VII. APPENDIX - LIST OF TRANSIENTS CONSIDERED FOR MATCHES 



Included below are transients with full spectral information. Transient XRFs are insuffi- 
ciently energetic for the current analysis. 

SWIFT GRB / Time / Trigger / RA / Dec / Duration T90 / Fluence (15-150 keV) / 1-sec 
Peak Photon Flux / Photon Index / XRT Initial Temporal Decay Index / XRT Spectral 
Index / (Gamma) Redshift / < Livetime > / Dl threshold 

050908 / 05:42:31 / 154112 / 20.451 (01:21:48) / -12.962 (-12:57:45) / 20 / 4.91 / 0.70 / 
1.86 (PL) / 1.33 / 3.9 / 3.340 / 0.92 / -0.389 

050730 / 19:58:23 / 148225 / 212.063 (14:08:15.1) / -3.740 (-03:44:24.0) / 155 / 24.2 / 
0.57 / 1.52 (PL) / 2.3 / 1.8 / 3.97 / 0.75 / -0.370 

050603 / 06:29:05 / 131560 / 39.982 (02:39:56) / -25.195 (-25:11:41) / 13 / 76.3 / 27.6 / 
1.17 (PL) / 1.78 / 1.71 / 2.821 / 0.92 / -0.385 

HETE GRB / Time Class / Coords / Redshift / Epeak / T90 / 30-400 keV fluence 
(erg/cm 2 ) / lightcurve / < Livetime > / Dl threshold 

GRB 020124A / 10:41:15 / XRR / 09h 32m 49.0s / -lid 27' 34" / 3.2 / 86.93 / 46.42 / 
6.1e-06 / 0.83 / -0.50 

GRB 020813A / 02:44:19 / GRB / 19h / 46m 37.9s / -19d 35' 16" / 1.25 / 142.1 / 89.29 
/ 8.4e-05 / 0.98 / -0.30 

GRB 020903A / 10:05:37 / XRF / 22h 49m 01.0s / -20d 55' 45" / 0.25 / 5 / 4.10 / 1.6e-08 
/ 0.92 / -0.22 

GRB 030429A / 10:42:22 / XRF / 12h 13m 06.0s / -20d 55' 59" / 2.66 / 35.04 / 10.30 / 
3.8e-07 / 0.95 / -0.352 

GRB 030723A / 06:28:17 / XRF / 21h 49m 29.7s / -27d 42' 07" / 2.3 / 11.3 / 6.78 / 
2.8e-07 / 0.96 / -0.352 
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TABLE V: GRBs recorded in southern hemisphere, with incomplete spectral or redshift information 



GRB 


Declination Angle 


< Livetime > 


Discriminator Threshold (V) 


051111A 


-2.64 


0.79 


-0.35 


051109B 


-53.959 


0.946 


-0.48 


051021A 


-45.519 


0.916 


-0.35 


051013A 


-4.399 


0.886 


-0.4 


051007A 


-0.512 


0.868 


-0.35 


051006B 


-67.715 


0.958 


-0.45 


051004A 


-10.869 


0.94 


-0.4 


051004B 


-15.347 


0.946 


-0.4 


051001A 


-30.4786 


0.928 


-0.47 


050927A 


-36.313 


0.934 


-0.44 


050921A 


-37.034 


0.928 


-0.43 


050917A 


-62.582 


0.97 


0.55 


0509 17B 


-21.242 


0.934 


-0.55 


050916A 


-50.5703 


0.976 


-0.57 


050915B 


-66.5899 


0.958 


-0.58 


050915A 


-27.9835 


0.97 


-0.58 


050911A 


-37.1506 


0.85 


-0.39 


050910A 


-29.775 


0.928 


-0.45 


050908A 


-11.0452 


0.838 


-0.39 


050907A 


-48.815 


0.898 


-0.39 


050830B 


-45.23 


0.694 


-0.36 


050826A 


-1.35683 


0.856 


-0.41 


050822A 


-45.9662 


0.736 


-0.41 


050818A 


-35.058 


0.97 


-0.51 


050817A 


-24.93 


0.916 


-0.44 


050813A 


-32.242 


0.718 


-0.4 


050813B 


-60.716 


0.412 


-0.4 


050801A 


-20.0719 


0.886 


-0.38 


050728A 


-26.984 


0.766 


-0.41 


050727A 


-21.768 


0.922 


-0.41 


050726A 


-31.9359 


0.952 


-0.41 


052726B 


-22.655 


0.904 


-0.41 


050725A 


-16.466 


0.898 


-0.42 


050724A 


-26.459 


0.904 


-0.40 


050721A 


-27.6189 


0.454 


-0.29 


050721B 


-1.44 


0.826 


-0.29 


050715B 


-0.059 


0.514 


-0.34 


050714A 


-15.528 


0.754 


-0.38 


050713A 


-12.264 


0.604 


-0.33 
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